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Abstract 

The natural outcome of theoretical calculations of microwave background 
anisotropy is the angular power spectrum as a function of multipole num- 
ber i. Experimental C^'s are needed for direct comparison. Estimation proce- 
dures using statistics linear in the pixel amplitudes as well as the conventional 
but less useful quadratic combinations are described. For most current ex- 
periments, a single broad-band power amplitude is all that one can get with 
accuracy. Results are given for the Capri-meeting detections. Mapping ex- 
periments, sensitive to many base-lines, can also give spectral "colour" infor- 
mation, either with a series of contiguous narrow-band powers or as parame- 
terized by a local "colour" index tiat (scale invariant is -2, white noise is 0). 
Bayesian analyses of the full first year DMR and FIRS maps give very similar 
band-powers {e.g., Qrms,PS = 17.9 ± c.f. 18.6ib4.7^K for toat = —2) 

and colour indices (with 1 and 2 sigma error bars) ti-at + 3 = 2.0^q4!^5 q '^^d 
l-8^o'g!l'j' 3 {c.f. the value 1.15 for a "standard" scale invariant CDM model). 
The 53 and 90 GHz DMR maps, as well as the FIRS map, have residual short- 
distance noise which steepens tt-at- This residual has so far been modelled 
by allowing the pixel error bars to increase, absorbing much of the effect, but 
further exploration is needed to see if a second residual evident in the data 
— which is, in part, responsible for the high hat — is from systematic errors 
or is physical. 
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"Un-cleaned" data with inadequate "dirt" models do not yet give us confidence that what 
is primary anisotropy has been well-separated from secondary backgrounds, foregrounds 
and instrumental systematics in any experiment. But we will. In the meantime, we need a 
phenomenology to display experimental results that (1) allows us to be unafraid of presenting 
data with the inevitable residual contaminants, (2) gives a AT/T estimator which is not 
sensitive to the specific experimental configuration (thus not rms anisotropics), (3) allows 
a meaningful comparison among experiments and (4) among theories. The angular AT/T 
power in broad bands that cover the multipole range explored by a given experiment satisfies 
these requirements. Whether the data allows one or many bands to be well-estimated 
depends upon the details of the experiment. 

Gaussian theories are completely characterized by the power spectrum, Cre = i{i + 
l){\aini\'^)/{2n), where ^(g) = E^m «<?m^£m(g) defines the aim- The signal {AT/T)p from 
the p^^ pixel of a CMB anisotropy experiment can be expressed in terms of linear filters 
^p,em acting on the a£m'- (AT/T)^ = J^im^p/mdem- The associated pixel-pixel correlation 
function can be expressed in terms of a quadratic Npi^ x Npi^ filter matrix Wpp// acting on 
Ct£, while its trace defines the average filter We {e.g., Bond 1989): 
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The "discrete-logarithmic integral" X of a function function is defined by eq. |^. Although 
there are as many as Npi^^Npi^ + l)/2 different filters to probe Ct£, in practice symmetries 



reduce this number to at most 
probes. 



Npix and often only a few are large enough to be effective 



1. Broad-Band Power and Angular Colour Measures 

If the power spectrum changes slowly over the band that = stretches over, one 
estimate of the broad-band power is just a renormalization of the experimental rms: 

{Ce)^ = iAT/T)l^JlWe] = I[WeCe]/I[We] . (4) 

For an experiment probing isolated (uncorrelated) pixels, this broad-band power is all you 
can get (since Cxpp' = C*Tii5pp')- Current intermediate and small angle experiments have in- 
sensitive Wpp'^e for pixel separations beyond the first few {i.e., Cxpp' approaches zero rapidly): 
splitting the M^-band into Bi,B2, ... sub-bands does little good. However, for mapping ex- 
periments like DMR and FIRS the number of pixel base-lines probed (number of useful 
Wpp',i) is large so the large H^-band can be split into a number of contiguous sub-bands 
(but not too many or the error bar matrix becomes too large and complicated). "Spectral 
colours" can be defined as logarithmic differences of the band-powers. Alternatively one 
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can assume a power spectrum with a local colour index tiat as well as a broad-band power 
amplitude (Q)^: 



Cm = {C^B + 1)'^""" ^Bi 1[W,]/1[W,{1 + . (5) 

Here an extra input estimator shape Ubi is included to give more flexibility for statistical 
testing. I usually choose Ubi = 1, although I find it is always instructive to see what a 
matched form, Ube oc Wi, gives. As we learn more experimentally about the true shape, 
Ubi can be matched to it. 

1 1 

The amplitudes {[(C£)^J 2 , [(C^)^_^] 2 ^ ...} can be determined by whatever statistical 
method we are most enamoured with, whether Bayesian or frequentist. Their likelihood 
function defines a multidimensional surface, whose peak gives the best band-power esti- 
mates and whose curvature matrix about the peak defines the error bar matrix; or, better, 
Bayesian credible regions can be defined for error estimation. 

It also turns out to be simple to transform old-style presentations of AT/T detections 
and limits to the band-power and colour language. 

The infamous (AT/T)c-curves as a function of coherence angle 6c can be translated 
into the more useful band-power estimates by noting the power spectrum for a "Gaussian 
correlation function" is of the form eq. ^ 

GCF : nAT = & Ube = e'^'/^ , u = {i + i)/(4 + i) , 4 + | = [2 sm{ej2)]-' ; (6) 
(Q)^ ^ {AT/T)l j[n2e-"'/2^,] /J[W,] . (7) 

Another form often adopted is the low-£ "Sachs- Wolfe" power for scalar metric pertur- 
bations {e.g., Bond and Efstathiou 1987), 



+ 1) (i-^y. (2 + 



It differs very little from the more natural form for phenomenology, eq. ^, with the index 
Hat related to the primordial density fluctuation index as shown. Thus = 3 is white 
noise in AT. 

For the DMR and FIRS beams (including pixelization effects), the relation between the 
quadrupole power C2 (or equivalently Qrms,ps) band-power is 

(^^)' - ^^2 - ^ {Ce)B e-"(--i)(^+°-^(--^)) , a,^r = 0.73 , aj.^s = 1.1 , (9) 

where Tq = 2.726 ±0.01 K. An advantage of the band-power over Qrms PS that the former 

is roughly independent of n^j-, while the latter is quite sensitive to it. Fig. 1 

Fig. 1 shows estimates for the power spectrum averaged over bands in £-space, assuming 
only wavelength- independent Gaussian anisotropics in AT/T are contributing to the (some- 
times cleaned) signals. See the reference list for the experiments' acronyms. In Fig 1(a), the 
data points denote the maximum likelihood values for the band-power and the error bars give 
the 16% and 84% Bayesian probability values (corresponding to ±la if the probability distri- 
butions were Gaussian). The horizontal location is at the average value {£) = T{iW)/2{W); 



3 



except for DMR, FIRS and Python, {€) is near the ly^-maximum. The horizontal error bars 
denote where the filters have fallen to e""'^ of the maximum. 

Proceeding from small £, Qdmr is the quadrupole power estimated using the 53GHz maps 
with a Galactic cut bccut = 20° (Bond 1993). The treatment of the dmr, firs and Tenerife 
broad-band powers are described later. The sp91 band-estimates are for the 9 point scan, 
the 13 point scan and for the combined statistical analysis of the 9 + 13 point scan (higher 
value). The offsets are for clarity. All 4 channels were simultaneously analyzed (Bond 1993). 
(The 4th channel of the 9 point scan gives a 1.5 x 10~^° 95% credible upper limit. Including a 
GCF with a synchrotron slope opens up the error bars of SP91 considerably.) The two recent 
MAX (M) results are for the scans in GUM (upper) and Mu Pegasus (lower, with a strong 
dust signal removed). The dotted hues ending in triangles denote the quoted 90% confidence 
interval for the MS AM single {g2) and double [gS) difference configurations for that half of 
their data that did not have obvious large sources in it. The dashed error bar shown next 
is the 84% upper credible limit for the SP89 9 point scan. The OVRO 7 point upper limit 
is last. For most experiments with detections, the 50% Bayesian probability (where the 
vertical and horizontal error bars cross) lies very close to the likelihood maximum. 

Fig. 1(b) gives band-powers for other experiments reported at the Capri meeting that I 
have estimated from GCFs (evaluated at or near the most sensitive angle 6c) or {AT/T)rms, 
but that I have not yet done a full Bayesian analysis for. The BigPlate result is denoted by 
bp, Python (py) is next, followed by Argo. These higher MSAM2 and MSAM3 results (than 
in Fig. 1(a)) are what one gets when the "source-half" of the data is also included. Finally 
the upper limit from the m = 2 mode analysis of the WhiteDish experiment is wd2. The 
dmr 2 number — with the smaller error bar than the dmr one in Fig. 1(a) — is a translation 
of the Bennett et al. (1994) result for the second year DMR data, Qrms,ps = 17.6 ± 1.7fj,K 
for = 1, to a band-power. Given my experience that transformed band-powers are in 
excellent agreement with the results of full Bayesian analysis, I am confident that Fig. 1(b) 
gives a fair account of the implications of the new detections reported at the Capri meeting. 

To show how well the transformation formulae Eqs. (^, |^, ^ work, I use some of the 
numbers quoted for the double difference Tenerife experiment (5.6° beam and 8.2° throw) 
at the Capri meeting. For the 15 and 33 GHz data, we heard that Qrms,ps = 26 ± 6^K for 
ris = 1, hence using (|) gives {Q)^ = 2.2tJ:^ x 10"^°. They also give' (AT), = 5At\t^K 
at Be = 4° (4 = 13.8). The ratio llu^e-'^^/We] /I[Wi] turns out to be 0.587, hence the 
transformation (|^) gives (Q)^ = 2.3^1'^ x 10~^°. For all the 15 GHz data (covering 700 deg^), 
Rebolo gave (AT/T)^^, = 39 ± lOfxK. I find I[We] = 0.98, hence {Ce) ^ = 2.11^9 x 10"^° 
using (^). Thus, all are quite consistent. 

Overlaying band-powers on theoretical Ci curves basically shows the status of the theories 
in question. However, since we do not know precisely where the observed power (Q)^^^^ 
lies in the band, the direct comparison should be with the theoretical band-power (Q)^^^ = 
X[CW]/2[W]. These are shown in Fig 1(c) for a few representative theories to illustrate the 
precision of power spectrum determination we need to finely differentiate among models. 
Error bars of 10%, the best DMR can possibly do, and that are not over-optimistic for 
intermediate- and small-angle mapping experiments, are also shown. To consider what is 
required for this accuracy, consider an experiment with Npi^ pixels with per-pixel error apix- 
Suppose the pixels are sufficiently separated that only Wi is an effective probe of For 
large Npi^, the z/-sigma uncertainty in the experimental value of the band-power is 
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To get 10% error bars requires experimental noise to be small c.f. the signal (we are basically 
already there), and Npi^^ = 200, i.e., a mapping experiment. For large Npi^, the observed 
maximum likelihood will fluctuate from (Q)^^^, the quantity we want, according to eq. [Tl|, 

but the error bars of eq. jl^ include these realization-to-realization fluctuations (thus a/2 
appears, not 1). 



2. The DMR and FIRS Colour Problem 

In this section, I describe some Bayesian results on the FIRS and flrst year DMR maps. 
The methods deal with all aspects of the maps simultaneously and so are highly sensitive to 
all components in it, whether they are the primary signals we are interested in or the warts 
that Dave Wilkinson admonishes us to be aware of and beware of, for they blemish "the 
face" that George Smoot (and we too) wish to contemplate. As we shall see, it seems to be 
pimples and not warts which I have to concentrate on cleaning up. Although I think even 
now I am doing this moderately well — for within the statistical observing strategies I use 
certain natural pimple-erasers are suggested — I still believe that the acne pattern colours 
my high precision look at the face. For the issue is one of colour. 

To make the DMR map data tractable for exploration, I used one lower resolution scale 
than the original maps, i.e., 5.2° rather than 2.6° pixels. The dipole and average subtractions 
were done after the Galactic latitude cut (taken here to be |6| > 25°), but before the 
lowering of the resolution. I used the Kneissl and Smoot (1993) revision of the DMR beam 
(which is non-Gaussian), and included corrections for digitization and pixelization, especially 
important in view of the 5.2° pixel size. Since COBE actually measures the difference 
between AT values at 2 beam-smeared points 60° apart, the pixel errors are correlated. 
A correction linear in the off-diagonal terms similar to one proposed by Lineweaver and 
Smoot (1993) has been applied to take this effect into account. The unknown amplitude in 
the average and quadrupole of the theory have been handled by treating them as marginal 
variables to be integrated over in the Bayesian analysis. 

The FIRS ('MIT') map (168 GHz) has a highly inhomogeneous weighting of each of the 
1.3° pixels. For the results shown here one lower resolution level was used, with 2.6° pixels, 
although I have checked that 1.3° pixels gives the same answer for the n^T = ~2 case. 
The average and dipole were removed in the map construction. To include the effects of 
pixelization and the best beam-smearing estimate, a 4.2° beam was used rather than 3.8°. 
(The differences in the amplitudes obtained are quite small.) 

A full Bayesian analysis of maps requires frequent inversion and determinant evaluations 
of Npix X Npix correlation matrices, the sum of all Crppi in the theoretical modelling plus 
the pixel-pixel observational error matrix C^pp'- I transform this matrix and the ampli- 
tude vector A = (Ap) to ones that absorb the pixel 'weighting', St = Cj^^'^CtCj^^'^ and 
X = Cjj^^'^A, with 'dimensions' of {signal — to — noisey'^, respectively. The removal of aver- 
ages, dipoles, gradients, etc. are included as more 'theoretical' signals, whose 'unobserved' 
contributions to Ap are integrated over, creating marginal likelihood functions. For their 
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prior distribution, very broad Gaussians arc adopted (essentially the same as uniform priors, 
but they regularize the inversions). I rotate to the eigenvectors of which are fully orthog- 
onal (uncorrclatcd) signal-to- noise modes for the map, with eigenvalues l + £TR,k and signals 
for each mode /c, where £tr = RStR^^ and ^ = Rx. With uniform weighting and all-sky 
coverage, these modes k are just the independent 5R(a£^) and ^'(a^m) for any theory, but 
with Galactic cuts followed by dipole removals they are complicated and theory-dependent. 
Even so, these modes are quite instructive linear combinations of the pixels. Often only 
a fraction of the modes are highly sensitive to theories being tested. For example, pixel 
differences within the beam are theory-insensitive, but highly susceptible to excess noise in 
the experiment, and the (5'/iV)-eigenmodes can be used for effective filtering or weighting 
strategies. 

I sort the (S'/A^)-eigenvalues, £TR,ki in decreasing values. These are plotted as open 
circles in Fig. 2 for DMR and FIRS, for the cases = 1 (a,h) and = 2 (c,d)\ the 1- 
sigma cosmic variance arc the two solid lines surrounding the points. £TR,k is normalized 
to (Tth — 1, where crf^j = (Q)^/10~^°; this just happens to be the amplitude that the data 
prefers. Although it may appear from the £TR,k decline that only the first few modes are 
effective at probing the theory in question, for uniform weighting and all-sky coverage it is 
basically equivalent to plotting Ci/{1{1+ 1)), which falls with ^ even without beam-smearing 
effects. A rough way to correct the visual impression is to multiply by k (which would be 
(£ -|- ly for homogeneous-weighting and all-sky coverage): k£TR,k is then nearly constant, 
decaying mostly because of the beam. 

The reason I didn't plot it this way is to make evident the (almost flat) excess power at 
high k in the data, which the beam-filtered theoretical points cannot possibly explain. The 
observational (solid) points in Fig. 2 are the average (5'/A^)-mode band-powers /?, 

^W^;J-E(S-1), ^/J^El- (12) 

'''P kei3 kep 

The (sparse)-binning was chosen to make the data trends clear. The error bars are estimated 
directly from the data, and include a correction for large scale power as well as a pixel noise 
error. 

The form of Fig. 2 was chosen to show a pictorial procedure for fitting the data: ath is 
adjusted until the error envelope on a^i^STR,k, which also scales with cr^, encompasses the 
data points as well as possible. Clearly it will not do very well for the high /c-bins for either 
the Hat = —2 or —1 theory. A first approximation to the residual 'noise', which is quite 
good for the high k end of the FIRS data, is to assume a constant offset. This can be viewed 
as increasing the error bars by a factor (1 -|- r) (so the offset is -|- 2r). Equivalently I 
can parameterize this residual noise by a new 'theory' source to be added to a^f^ £TR,k-i with 
(5'/A^)-power (tI^sStr = r'^ + 2r, where £tr = N~lJ2k^TR,k, arranged so that ff ares = crth, 
the power would be equal. 

Of course, this pictorial procedure with quadratic combinations of the modes is not the 
best statistical approach to deal with the data. Once one has transformed to the (S/N)- 
eigenmode system, a Bayesian analysis on the full map is very straightforward and fast, even 
with average and gradient subtractions. The effect of these subtractions is hard to show in 
Fig. 2: the /c-modes are no longer statistically independent. (The average affects the first 
bin, but the dipole has influence out to /c ~ 30, and is responsible for some of the power loss. 
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Note that the observational /9-powers can become negative.) The first step in the Bayesian 
method is the construction of a joint likelihood function in ath and (Tres- The contour maps 
in Fig. 3 show a strong maximum in ath at almost the same level for the DMR 53A+B, 
90A+B, and FIRS maps, in spite of the large differences in the position of the maximum in 

The single broad-band powers for DMR and FIRS plotted in Figure 1 arc derived from 
the marginal distributions in signal amplitude, found after integrating the joint distribution 
over all possible residual noise amplitudes {i.e., integrating over the of Fig. 3). 

For Us = 1, Qrms,ps = '^T-QfJ-K {{Ci) ^/10~^^y. I shall express the results in the Qrms,ps 
terms: for the Ssi + B GHz map, I get 17.9 ± 2.%ixK (17.6 ± 2.8^^-?^ if I ignore off-diagonal 
terms in Copp', 17.2 ± 2.5fiK if I use a 7° Gaussian beam with no pixclization correction). 
This compares with the DMR team's best estimation for the first year data of 17.1 ±2.9fiK 
(Wright et al. 1994) and the new two-year result of 17.6 ± 1.7 fiK. For 90A-I-B, I get 
18.6 ± S.SfxK and for 31A+B, 15.8 ± 5.8fxK. For 53A-B, I get Qrms,ps = ^.StUfiK, while 
for 90A-B, O.OtoioA*-?^: i-e., no spurious large scale power in the difference maps. (Guriously 
I find no clear detection in 90A; however the error bars for this map are quite large.) For 
Us = 1 and the FIRS map, the level is IS.Gll^^^ixK - i.e., at the DMR level. 

The probability for the residual noise level is found by integrating over the crt/i-axis. 
The most probable value is indicated by the horizontal dotted line joining the triangular 
points in Fig. 2, about at the level the eye would pick. For FIRS, zero in residual noise is 
excluded at more than the 10 sigma level (with r = 0.25 most probable, independent of tzat 
(Fig. 2(c) c.f. Fig. 2(a)). However, the high plateau around k = 200 in Fig. 2(a,c) suggests 
0"pj2:-enhancement is not the whole story (in accord with failure by the FIRS team to identify 
the source of the residual after exhaustive checks). 

For the first year DMR data, the residual depends upon the map and upon the spectral 
steepness. I find no evidence for a residual in the 31A-I-B map, but do have them in both 
53A+B and 90A+B. The residuals in 53B (r = 0.07) and 90B (0.09) are larger than in 53A 
(0.01) and 90A (0.02). 53(A-B) has residuals about the same as 53A+B, suggesting that it 
is not the effect of physical sources on the sky. The numbers in brackets are for uat = —2. 
They drop with increasing uat (see Fig. 2 for tiat = ~1)- 

Al Kogut and the DMR team have checked the scatter in their data before conversion 
to a map and confirm the trends I find in r, but with lower amplitudes than my n^T = 
—2 values, more like my uat = ^1 values, the index the data prefers. However, there 
are effects not usually included in map analysis which will increase the variance, e.g., the 
"unobserved" Galactic plane pushes up the variance in off-plane pixels through the 60° 
correlation. Modelling the excess in Fig. 2(c,d) by a constant is clearly not optimal. The 
shape of the data points is suggestive of beam-effects not being properly included. However, 
even with the old standard 7° beam (quite different than the Kneissl and Smoot 1993 beam), 
and with no pixelization corrections, the SxR^k curve still falls off faster than the data. 

Spectral indices from uat = — 3 to uat = in steps of 0.1 were run to construct 
likelihood functions in {^at, crt/i, CTres}- Integrating first over ares (marginalizing it) allows 
one to construct nAT-<^th contour maps, not shown here. Integrating again, over ath^ gives 
the probability distribution for uat- The colour index I get is uat + 3 = 2.0lo;4;l'j'^o f*^^ ^^e 
DMR53A-I-B map and 1.8l2:^:t?| for the FIRS map, with the second errors given denoting 
2-sigma values. The reason for such steep tiat is evident from Fig. 2: higher uat gives a 
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slower fall of the mode power with increasing k. However, one cannot increase n^r too much 
or the flattening of SxR^k at low k becomes too much for the high amplitude data points 
there to contend with: uat = — 1 is the compromise. 

It might be thought that filtering out the high k modes will allow us to avoid modelhng 
the excess found there. However, so far I find that sharp /c-filtering, keeping only a few 
hundred modes, does not help for FIRS. The k ~ 200 plateau keeps the preferred index 
high. For DMR with k < (16)^, the residual r grows for low uat, stealing some of the 
power from ath (which drops by 20% for uat = —2). However, the n^T value (and error 
bars) remain about the same as when all of the modes are included. For k < (13)^, uat 
becomes poorly determined, l.OlJ^y^^^o- ^ faulty beam structure is responsible, adding to 
the likelihood function a beam-smearing parameter could be used to set it by finding the 
most probable value. However, a 7° Gaussian beam and no pixelization corrections gives 
1.7lo!4^o;9, still steep c.f. the inflation prediction. Models with three or more components 
have yet to be explored. A particularly interesting case is a white noise uat — distribution 
to model small scale (non-primary?) sources, plus a large scale power source with variable 
n-AT, plus pixel error bar enhancement. 

Most people probe the data using quadratic combinations Qa of the pixel amplitudes 
(AT/T)p, and hence of ^k- Qa = J2kk' Ikk'^kCk' ■ The correlation function and various 
combinations for the C^-spectrum are examples, as, of course, are the (5'/A^)-powers Vs/n,i3- 
The sum of the Vs/N,k hnearly weighted by the (5'/A'")-eigenvalues, Y.k ^TR,k^k^ Boughn- 
Cottingham (BC) statistic. Its effect can be seen in Fig. 2: multiplying the data by STR,k 
strongly suppresses the contribution of the deviant high fc-modes to the sum. A pictorial 
exercise such as in Fig. 2, but with only one data point, can be done for the BC statistic. 
A better method is to construct a likelihood using many (signal-plus-noise) Monte Carlo 
realizations of the ^k, with average and dipole subtraction. This is the map. Thus, map 
construction is very fast since it involves just Npi^ random number choices and no direct 
summation over Yim at each pixel. This high speed helps, since to create a reasonably 
smooth likelihood function (rather than just using frequcntist statistical measures) requires 
fifty thousand rather than a few thousand realizations. The eigensystem is of course also 
excellent for sets of quadratic statistics, e.g., for angular bins of the correlation function and 
^-bins of various prescriptions for quadratic ^-space power spectrum estimators. One just 
has to rotate the pixel-pixel quadratic operators into the 7^^, form in A:-space. 

The BC quadratic allows good recovery of the most probable Bayesian values (but it costs 
more in computer time since no Monte Carlo runs are required for the Bayesian calculation). 
With appropriate normalization, the BC statistic turns out to be the most sensitive single 
quadratic estimator of the band-power if there were no dipole etc. removals and no other 
signals in the data. For FIRS, the BC value is a little bit high. I have therefore constructed 
best single band-power estimators with these effects included. For FIRS, the (Ci)^ result 
is bang-on the Bayesian one. For DMR53A-I-B, both the BC statistic and my modification 
give the same value, lower than the Bayesian value with all modes included and closer to 
the k-cut value described above: the two forms of filtering are not that different. (Wright et 
al. 1994 have also used the BC quadratic on the DMR data, and find similar amplitudes.) 
Unfortunately, while well-suited to Qrms,ps estimation, these quadratic forms change with 
uat-i so joint likelihoods should not be constructed, and, in any case, are a poor way to get 
at Uat- The obvious first thing to try is Weiner-filtering (with £TR,k/i)- + £TR,k) acting on 



8 



^fe); however, not only does it change form with n^T, it also changes form with ath, and 
so it is not useful for either Qrms,ps or "at statistical analysis, although it can be used to 
construct cleaned maps. I have also tried Bayesian analyses without a residual, but with 
the mode sum in the log-likelihoods suppressed by a measure oc — 1 gives results in 

agreement with the BC-style quadratic statistics, but with smaller error bars because it deals 
with the entire map. Unfortunately, the result is not stable to variations in a, especially for 
FIRS, since lower a does not sufficiently suppress the residual. 

It is impractical to construct a likelihood when there are a number of 7"^'s to compare 
with: to make smooth functions, too many Monte Carlo runs are required even for this 
diagonal frame. Usually some distance measure of observed from simulated Qa is used {e.g., 
the simple chi^ measure for C{9) used by the DMR team in Smoot et al. (1992), and variants 
which include the effect of the strong correlation from angle-bin to angle-bin in the theory, 
which I have used Bond 1993). 

What the do is filter pixel-pairs, and this can have great advantages over just linear 
filtering of modes. For example, although I do not know the details of the angular structure 
of the residual noise, the pixel-error enhancement model contributes only to the zero angle 
bin of the correlation function. If the residual's effect spills over to neighbouring pixels, one 
can just cut out the core of the correlation function, say out to Of^hm- This would be an 
excellent filter if some part of the residual arises because the beam is not quite right or if 
there are white-noise or Poisson sources on the sky. The same effective filtering does not 
happen with quadratic band-power estimators. For example, white sky-noise gives a rising 
£^ contribution to C^. The residual in the FIRS map obscured the interpretation of my 
results using quadratic band-power estimators (which differ from the Hauser-Peebles form 
used by Wright (1993) for DMR, which might be better). Quadratic band-power operators 
which exclude pixel-pairs out to 9fu,hm do filter this rise, but the interpretation as power in 
an £-band is lost. Thus, I believe that the Bayesian approach is superior to using a quadratic 
set of band-power operators for C^-estimation. And what comes out is an estimate of the 
ensemble-averaged power, which is what we want, and not the specific realization of the 
power that exists in a map, even if that map is our own sky. 

The filtering of small-angle systematic or physical effects in correlation function analysis 
may make the n^T determined with C{9) a better indication of the angular colour of the 
large-angle sky than I get - at least until the nature of the high A;-power is better understood. 
For the FIRS map, there are indications from C{6) for a shallower slope, uat ~ —1.7 
(Ganga and Page 1994); even so, it is clear from Fig. 2 that there are some anomalies in 
the S/N band-powers that are not favourable to uat < —2. Bennett et al. (1994) find 
n^T + 3 = 1.59lg:i for two years of DMR data. Modifications of the Bayesian approach to 
model a more extended residual noise structure as well as the single-pixel level one is the 
next step for me to better observe the angular colour of the CMB sky. 

I would like to thank the FIRS and DMR teams for much discussion on how to treat their 
maps, in particular Lyman Page, Al Kogut and George Smoot. This research was supported 
by NSERC at Toronto, CNRS in Paris, and the Canadian Institute for Advanced Research. 
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FIGURE CAPTIONS 



Figure 1. (a) and (b) give band-power estimates for the vintage Fall 1993 data. The labels 
for the experiments are defined in the text. The vertical error bars are 1 sigma, the inverted 
triangles are 90 or 95% confidence limits. Systematic errors or non-primary sources {e.g., 
dust, synchrotron) may contribute to these C^'s; the underlying primary Ce may be lower, 
but could be higher because of 'destructive interference' among component signals. Ci is 
shown for a "standard" CDM model (f2 = 1, = 1, h = 0.5, Qb = 0.05) with normal 
recombination to illustrate a typical Doppler peak, but the early reionization (at z > 200) 
Ce is decidedly peak-less. Whether a Doppler peak exists is unclear from the current data, 
(c) Theoretical band-powers with (eventually) achievable 10% error bars are displayed for 
both the Hs = 1 CDM model and a chaotic-inflation inspired Us = 0.95 one - including a 
gravity wave contribution (Crittenden et al. 1993). The dotted Ci has 7^ and high Hq 
(i^A = 0.75, Vtcdm = 0.22,^5 = 0.03, ifo = 75, = 1). This shows the high precision we 
require to differentiate among models (for more details, see Bond et al. 1994a,b). 
Figure 2. (S/N) band-powers Vs/n,i3 (eq. for the orthogonal signal-to-noise eigenmodes 
of the FIRS and DMR 53A+B maps, and for uat = —2, —1 powers, as indicated. The error 
bars denote observational variance for the quadratics. The small open points denote the 
/?-band averages of £TR,k- The upper and lower curves about these give the theoretical 
variances. All three curves scale with cr^^^. The triangles joined by the dotted curve show 
the most likely residual noise the Bayesian analysis finds. The observed Vs/n,i3 should be 
matched within the errors by a^f^STR,k + crles^TR through adjustment of ath and the residual 
ares- For the DMR cases, the dashed line shows Etr^u for a 7° Gaussian beam with no 
pixelization correction. The data points also move a bit, but not enough for this beam to 
explain the high fc-behaviour. Small k are, of course, unaffected. Although cr^es drops, ath 
hardly does. 

Figure 3. Contour maps (1-sigma to 10-sigma) of DMR and FIRS likelihood functions 
for the scale invariant tt-at + 3 = 1 case. The heavy dot denotes the maximum. Here, 
o"t/i = [{C't) bI'^^^^^Y^'^ ^"^^ residual noise amplitude ares is normalized to have the same 
total power as in the CTf/^-band if ares = crt/i- The 53A+B, 90A+B and FIRS maps show 
the same band-power detection, within the errors, and in spite of the differing residuals. 
The 53 GHz difference map, 53A-B, shows no large scale power, although the residual offset 
remains at the same level as in A-l-B. The same holds for 90A-B. As tzat increases, the 
required residual goes down for the 53A+B map, ares ~ 0.466v^— nAr, with no need of a 
residual for white noise (jt-at = 0), whereas it remains constant at ares ~ 2.5 for the highly 
inhomogeneously-sampled FIRS map, . 
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